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BUCKLING OF ECCENTRICALLY STIFFENED ORTHOTROPIC 


CYLINDERS UNDER PURE BENDING 

By David L. Block 
Langley Research Center 

SUMMARY 

The stability of stiffened orthotropic cylinders which include the effects of eccentric 
stiffeners and are loaded with pure bending or any combination of bending and compres- 
sion is investigated analytically by means of a small deflection theory. Solutions to the 
buckling equations are obtained for simple support boundary conditions by use of the 
Galerkin method. The pure bending buckling loads are compared with existing pure com- 
pression buckling loads for three contemporary cylinder configurations. Also shown are 
predicted buckling modes and typical interaction curves between bending and compression. 
The results show that eccentricity effects are substantial and that the maximum pure 
bending load may be as much as 40 percent greater than the pure compression load. 

INTRODUCTION 

The effect of eccentricity or one-sidedness of stiffening elements in determining 
buckling loads of stiffened cylinders has been shown to be of great importance. (See 
refs. 1 to 4.) In reference 3, a theory for eccentrically stiffened orthotropic cylinders 
employing Donnell-type assumptions and neglecting prebuckling deformations was derived 
from energy principles, and buckling solutions were obtained for a simply supported cyl- 
inder under any combination of uniform axial load and lateral pressure. 

The purpose of this paper is to extend the analysis of reference 3 to include the 
buckling of an eccentrically stiffened orthotropic cylinder subject to a loading of pure 
bending or combination bending and compression (nonuniform axial loading). The equa- 
tions developed are also valid for an unstiffened cylindrical orthotropic shell. Solutions 
to the buckling equations corresponding to classical simple support boundary conditions 
are obtained by use of the Galerkin method in a manner similar to that employed in refer- 
ence 5 in which an unstiffened isotropic cylinder under pure bending was investigated. 
Maximum compressive stress resultant at buckling due to bending, hereafter referred to 
as maximum bending loads, are compared with uniform compression loads for three 
types of cylinders with contemporary geometry, that is, ring- stiffened corrugated cylin- 
ders, ring- and stringer-stiffened isotropic cylinders, and longitudinally stiffened 



isotropic cylinders. Interaction curves between bending and compression and the buck- 
ling mode shapes are also presented. 


SYMBOLS 


The units used for the physical quantities in this report are given both in the 
U.S. Customary Units and in the International System of Units (SI). The relationship 
between these two systems of units can be found in reference 6. 

A cross-sectional area of stiffener 


a n>^n ,c n 


D x ,Dy 


D xy 

d 

E 

E X ,Ey 


length of stiffened cylinder (see fig. 1) 
constant coefficients of deflection functions 

bending stiffnesses of orthotropic plate in longitudinal and circumferential 
directions, respectively 

twisting stiffness of orthotropic plate 
stringer spacing (see fig. 1) 

Young's modulus 

extensional stiffnesses of orthotropic plate in longitudinal and circumferential 
directions, respectively 

shear modulus 


J xy 


x o 

i,j,n,q,p 


J 

l 


m 


N„ 



inplane shear stiffness of orthotropic plate 

moment of inertia of stiffener about middle surface of shell 

integers 

torsional constant for stiffener 
ring spacing (see fig. 1) 

number of half waves in cylinder buckle pattern in longitudinal direction 

total compressive stress resultant 

maximum compressive stress resultant due to bending 

critical maximum compressive stress resultant due to bending when N Xc = 0 
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uniform compressive stress resultant 



N v critical compressive stress resultant when N Y =0 

x c x b 

R radius of cylinder to middle surface of shell (see fig. 1) 


t 

l 

u,v,w 


thickness of isotropic cylindrical shell wall 

A 

effective wall thickness of stiffened isotropic cylinder, — S- + t 

d 

displacements in x-, y-, and z-directions, respectively, of a point in middle 
surface of shell 


x,y,z orthogonal curvilinear coordinates with origin lying in middle surface of shell 

(see fig. 1) 

z distance from centroid of stiffener to middle surface of shell (see fig. 1), 

positive if stiffener lies on external surface of shell and negative if stiff- 
ener lies on internal surface 

Kronecker delta, 6^ = 0 when i * j and 6^ = 1 when i = j 

p v ,ju Poisson's ratios for bending of orthotropic plate in longitudinal and circum- 

x y 

ferential directions, respectively 

/i v ' ,p ' Poisson's ratios for extension of orthotropic plate in longitudinal and circum- 

x y 

ferential directions, respectively 


Subscripts: 


s,r stringers (longitudinal stiffening, parallel to x-axis) and rings (transverse 

stiffening, parallel to y-axis), respectively 

x,y longitudinal and circumferential directions, respectively 

A subscript preceded by a comma denotes partial differentiation with respect to the 
subscript. 


THEORETICAL SOLUTION 

The equations of equilibrium governing the buckling of a stiffened cylinder subject 
to inplane axial stresses only are (ref. 3): 


E„ 


1 - JU 'JU ' 

^x *y 


u ’xx + 


^s^s 


( u, xx ” z s w, xxx 


+ 1 - 1*1 ' ( v, *y + it) + G *y( u ’yy + v ’xy) = 0 

x y (la) 
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i - r I i \’yy r 1 ’yyy/ 

x y 


Mx'E 


+ 1 - m x ' M y ’ u,xy + Gxy ( u ’ xy + v,xx ) 0 


(ib) 


D x EsI °s 

(1 nr^T + “d^! W ’xxxx + \l . ^ ■ xy ■ 1 - Vy 


MyE x OT ~ Mx^y G S J S . G r J r 

+ 2D + . + d l i ■■ Jxxyy 


w >. 


D„ 


E r I 


r*o T 


w 

\i - uu ■ * ] w, yyyy + R /i . n A^ x u ’ x + v, y + R 


x' y 


r(i - m x * My ') 


E S A 


E r A r _ (2w,, 

d Z S U >XXX - Z v\ R 


S^S ; 




When equations (1) for a stiffened cylinder (fig. 1) were developed, the following 
basic assumptions were made. The stiffened cylinder is composed of an orthotropic shell 
with infinite transverse shear stiffnesses and with elastic constants as given in refer- 
ence 7. The stiffening elements are uniform equally spaced rings and stringers which 
are assumed to have elastic properties and to be closely spaced so that their elastic prop- 
erties can be averaged over the stiffener spacing. The usual Donnell-type assumptions 
are used to specify buckling displacements in the shell. The stiffeners are treated as 
beam elements with stiffener twisting accounted for in an approximate manner. In cases 
where both rings and stringers lie on the same surface of the shell, the effect of joints in 
the stiffening framework is ignored. 

The cylinder loading is taken to be a combination of bending and compression, in 
which the inplane loading prior to buckling is given by 

N =N +N x cos R (2) 

x x c b K 

where N x is the uniform compressive stress resultant and N x ^ is the maximum 
bending stress resultant; both are positive in compression. 

The boundary conditions to be satisfied are those of simple support which are given 
in reference 3 as 

w = 0 
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D, 


Eel 


1 " M x M y 


/ % ‘ t ' sA °S 

( w -xx + v-yy) + -T^ w ’: 


EoAc 

, - z u, = 0 

XX d S ’X 


E, 


1 - u 'u ' 
^ x ^ y 


u ’x + ^y\’y + 1 


E C A C E q A c 

+ - -VV>xx = 0 


v = 0 


Expressions for the displacements which satisfy these boundary conditions are: 


u = cos 


nijrx 


I «n«-f 


n=0 


v = sin2 l7rX 


Y b n sin SSL 

L n R 

n=0 


w = Bln2!Ii £ a n cosffi 


n=0 


(3) 


where m denotes the number of axial half waves in the buckle pattern. 

By substituting equations (3) into equations (la) and (lb), the coefficients b n and 
c n may be determined as functions of the coefficients a n and from equations (lc), (2), 
and (3) the following equation is obtained: 


oo 

Q(w) = y a n 


n=0 


a . A ,1^13 - A 11 A 23\ . A l A 12 A 23 " A 13 A 22\ 

Ao ° + A 23l-r^ r~r + A m~-— — — 2 - 


33 


A 11 A 22 " A 12 


A 11 A 22 _ A 12 


(> 


N x + N x cos Ife) 2 
x c x b rA a 


sin 2^ cos ^-=0 


in which 


A 


11 “ 


l 1 - M X V y’ 


E S A S1 




* 


(4) 
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where the reciprocal relations m *E = fx 'E and ju. D = u D (ref. 7) have been 

xy yx x y yx 

employed to simplify the equation. 

By use of the Galerkin method, a set of equations for the unknown coefficients a n 
is then 

r% 3. r\ 2ttR 

\ \ Q(w) sin E— ^ cos ^ dx dy = 0 (5) 

J 0 Jq a R 


from which a homogeneous set of equations is determined as 


a J F - N x 
n\ n x c 




In " u 0n) a n- 1 + a n+l 


1 + a- - 6, 


3 - 


( 6 ) 


where 


F = (— Y 
n \mir/ 


A A , A 1 2 A 1 3 ‘ A 11 A 23\ a / A 12 A 23 “ A 13 A 22\ 
a 33 + A 23 ~ : : 1 ~) + a 13 


A 11 A 22 " A 12^ 


L 11 A 22 " A 12 
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n = 0, 1, 2, 3, . . 


V 1 

(j = n) 

6 jn = ° 

(j *n) 


Equation (6) can be used to determine buckling loads under pure bending and pure 
compression (where the equation reduces to the result given in ref. 3) and to investigate 
interactions between bending and compression. For nontrivial solutions of equation (6), 
the determinant of the coefficients must equal zero. When computing the buckling load 
for a specified loading, the determinant is minimized numerically for integral values 
of m, the number of axial half waves in the buckle pattern. The size of the determinant 
is determined by requiring that the buckling load converges to a desired accuracy. For 
the calculations presented herein, the size of the determinant was increased until the 
buckling load between a specified order determinant and a determinant of one less order 
did not vary by more than 0.5 percent. All the numerical calculations presented herein 
were performed on an electronic computer. 

RESULTS AND DISCUSSION 

Because of the large number of geometric parameters involved in equation (6), it is 
impractical to present results of a general nature. However, it is of value to present 
some computed results for cylinders of contemporary proportions in order to study 
loadings of pure bending compared with pure compression, eccentricity effects under 
these loadings, and interaction between bending and compression. Therefore, computa- 
tions of bending loads were made for three types of stiffened cylinders appropriate for 
large-diameter booster interstage structures: ring -stiffened corrugated cylinders, ring- 
and stringer-stiffened isotropic cylinders, and longitudinally (stringer) stiffened isotropic 
cylinders. For the calculations, the material in the cylindrical shell and the material in 
the stiffeners were taken to be identical, and a value of 0.32 was assigned to Poisson's 
ratio for the stiffened isotropic shells. The dimensions of the cylinders are given in fig- 
ure 2. Eccentricity effects for each cylinder configuration were studied by moving the 
stiffeners from the internal to the external surface of the shell. 

The results of the computations are presented in table I and in figures 3 to 9. 
Details and discussion of the computed results for each cylinder configuration follow. 

The maximum pure bending loads for general instability failure of ring-stiffened 
corrugated cylinders and for general and panel instability failure of ring- and stringer- 
stiffened cylinders are presented in table I and are shown in figures 3, 4, and 5, respec- 
tively. In the figures, a nondimensional buckling load parameter N X] JE X or N x j E x 
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for corrugated cylinders and N Xb /Et or N^Et for ring- and stringer-stiffened 
cylinders, is plotted against a nondimensional ring spacing l/ R. For the panel insta- 
bility calculations, the ring spacing l was taken to be the length of the cylinder, and 
thus the cylinder is a longitudinally (stringer) stiffened isotropic cylinder. Results are 
presented for cylinders loaded in pure bending (computed by using eq. (6)) or loaded in 
pure compression (computed by using eq. (15) or (17) of ref. 3). 

For the ring- stiffened corrugated cylinders studied, the computed results (table 1(a) 
and fig. 3) indicate that eccentricity effects are as important for pure bending loads as 
those observed for pure compression (ref. 3). For internally stiffened cylinders, the 
critical maximum bending load was found to be as much as 40 percent greater than the 
critical pure compression load. (See last column of table 1(a).) This percentage increase 
is somewhat greater than similar results of 25 percent presented in reference 4; how- 
ever, the corrugated cylinders in reference 4 do not have the same properties as the ones 
examined herein and, therefore, the two results would not be expected to agree. 

For the ring- and stringer-stiffened and the longitudinally stiffened isotropic cylin- 
ders, the computed results (tables 1(b) and 1(c) and figs. 4 and 5) show that the differences 
between external and internal stiffeners are not as large as those observed for the corru- 
gated cylinders, but nevertheless, both effects may be substantial and must be considered. 

In figures 6 and 7, typical interaction curves between bending and compression are 
shown for the corrugated cylinder and for the ring- and stringer-stiffened cylinder, 
respectively. The interaction curve for the internally stiffened corrugated cylinder pos- 
sesses a vertical tangent at the abscissa as observed in reference 4; however, for the 
other interaction curves this behavior was not as pronounced. Hence, except for the 
internally stiffened corrugated cylinder, the interaction curves presented may be approx- 
imated by straight lines. 

In figures 8 and 9, the buckling mode shapes for pure bending and pure compression 
are shown (corresponding to the end points of the interaction curves of figs. 6 and 7). The 
buckle shapes for pure bending show that the buckles are localized on the compression 
area of the cylinder, a result observed in references 4 and 5. 

In addition to these cylinder configurations, calculations were made for unstiffened 
isotropic and orthotropic cylinders by setting the ring and stringer properties equal to 
zero and using the appropriate cylindrical shell constants in equation (6). For isotropic 
cylinders, the calculations showed that the effects of bending are negligible when compared 
with compression which agreed with references 5 and 8. The three cylinder configurations 
considered herein were analyzed as orthotropic shells by using orthotropic constants 
defined in reference 3. The calculations for these cylinders showed that even though 
eccentricity effects are neglected, there can be a considerable increase in the critical 
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compression load due to bending; for example, the corrugated cylinder gave an increase 
of 18 percent. 

Additional information can be obtained by considering the circumferential buckling 
modes for pure compression which are given in reference 3. The results presented in 
table I indicate generally that the differences between pure bending and pure compression 
loads may be considered small when the cylinder has an axisymmetric pure compressive 
buckling mode (for example, externally stiffened corrugated cylinders or unstiffened 
isotropic cylinders) or when the cylinder has a pure compressive buckling mode with a 
large number of circumferential waves (longitudinally stiffened cylinders). However, the 
results indicate large differences between bending and compression for cylinders which 
have pure compressive modes consisting of intermediate numbers (for example, 6 to 15) 
of circumferential waves, such as the internally stiffened corrugated cylinders or the 
ring- and stringer-stiffened isotropic cylinders. 

Also from table I it may be noted that the size of the determinant needed for con- 
vergence will depend upon the number of waves in the circumferential pure compressive 
buckling mode; that is, cylinders buckling with a small number of circumferential waves 
in compression will require fewer terms to approximate the buckle shape in bending and 
a smaller determinant order will be required for convergence of the buckling load. 

CONCLUDING REMARKS 

A stability equation for eccentrically stiffened orthotropic cylinders subject to a 
loading of pure bending or to combination bending and compression has been derived by 
means of the Galerkin method for simple support boundary conditions. Sample calcula- 
tions have been shown for three types of contemporary stiffened cylinders to illustrate the 
importance of bending loads, eccentricity effects under these loads, and interaction 
between bending and compression. The calculations presented show that the difference 
in predicted buckling load between a cylinder loaded in pure bending or pure compression 
may be considered small when the pure compressive buckling mode is axisymmetric or 
has a large number of circumferential waves. The buckling analysis of cylinders, which 
buckle in compressive buckling modes with a moderate number of circumferential waves 
and which are to be loaded in bending, require the use of an analysis which accounts for 
the bending load. The calculations also show substantial eccentricity effects in bending 
comparable with those obtained for cylinders in compression. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 14, 1966. 
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TABLE I.- BUCKLING CALCULATIONS FOR STIFFENED CYLINDERS 


(a) Ring-stiffened corrugated cylinders (b) Ring- and stringer -stiffened cylinders (c) Longitudinally stiffened cylinders 


1/ R 

N Y /El 
x b/ 

Minimum 

determinant* 

N v /N v 
x b / x c 

Stringers, e 

xternal; 

rings, external 

0.05 

0.006239 

5 

10 

1.045 

.10 

.005318 

4 

12 

1.038 

.15 

.004972 

4 

12 

1.050 

.20 

.004788 

4 

12 

1.072 

.25 

.004666 

3 

14 _____ 

1.089 

Stringers, e 

xternal; 

rings, internal 

0.05 

0.004752 

4 

14 

1.156 

.10 

.004335 

3 

14 

1.130 

.15 

.004144 

3 

14 

1.114 

.20 

.004025 

3 

16 

1.109 

.25 

.003938 

3 

16 

1.101 

Stringers, i 

nternal; rings, internal 

LO 

o 

o 

0.004470 

3 

14 

1.187 

.10 

.003949 

3 

14 

1.151 

.15 

.003669 

3 

14 

1.159 

.20 

.003476 

3 

16 

1.163 

.25 

.003319 

2 

16 

1.194 


l/R 

%/ E x 

m 

Minimum 
determinant * 

N x /n t 

x b / x c 

Rings, external 

0.15 

0.005514 

2 

6 

1.089 

.20 

.004679 

2 

6 

1.091 

.25 

.004171 

2 

6 

1.090 

.30 

.003828 

2 

6 

1.088 

.35 

.003581 

2 

6 

1.086 

.40 

.003391 

2 

7 

1.083 

.45 

.003245 

2 

7 

1.080 

.50 

.003128 

2 

7 

1.079 

Rings, internal 

0.15 

00 

ID 

a 

o 

© 

o* 

1 

10 

1.396 

.20 

.002084 

1 

10 

1.397 

.25 

.001839 

1 

11 

1.389 

.30 

.001666 

1 

11 

1.376 

.35 

.001536 

1 

11 

1.376 

.40 

.001442 

1 

11 

1.352 

.45 

.001353 

1 

11 

1.341 

.50 

.001286 

1 

11 

1.329 


l/R 

N * b /Ei 

m 

Minimum 

determinant* 

N x b /N Xc 

Stiffeners, external 

0.15 

0.007133 

1 

12 

1.013 

.20 

.004810 

1 

12 

1.015 

.25 

.003681 

1 

26 

1.019 

.30 

.002866 

1 

32 

1.018 

.35 

.002291 

1 

34 

1.022 

.40 

.001883 

1 

36 

1.022 

.45 

.001596 

1 

34 

1.027 

.50 

.001379 

1 

34 

1.029 

Stiffeners, internal 

0.15 

0.006888 

1 

30 

1.021 

.20 

.003909 

1 

28 

1.027 

.25 

.002522 

1 

28 

1.032 

.30 

.001769 

1 

28 

1.037 

.35 

.001315 

1 

28 

1.047 

.40 

.001022 

1 

26 

1.048 

.45 

.000821 

1 

26 

1.053 

.50 

.000679 

1 

26 

1.058 


^Minimum order of determinants required for 0.5-percent error in buckling load. 
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R = 200 (508) 


p =4.00(10.16) 

9 = 600 

t * 0.16 (0.41) 
c 


0.10 

(0.25) 



(a) Ring-stiffened corrugated cylinders; a = 200 in. (508 cm). 


0.10 

(0.25) 



R - 200 (508) 


Stringer 


0.90 _J 
(2.29) 1 

— 

0.10 „ 

t 

2.00 

(0.25) 

(5.08) 

!— i l 

L_ 0.90 
1 (2.29) 


Stringer section 


(b) Ring- and stringer-stiffened cylinders; a = 200 in. (508 cm). (See fig. 2(a) for ring section.) 

Figure 2.- Dimensions of stiffened cylinder walls. Dimensions are in inches. (Parenthetical dimensions in cm.) 
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Maximum pure bending 



Figure 3.- General instability predictions for ring-stiffened corrugated cylinders. 





Figure 4.- General instability predictions for ring- and stringer-stiffened cylinders, 







Pure compression (ref. 3) 
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Figure 9.- Predicted mode shapes for ring- and stringer-stiffened cylinder with i/R - 0.15, 



Figure 4.- General instability predictions for ring- and stringer-stiffened cylinders, 



Maximum pure bending 



Figure 5.- Panel instability predictions tor longitudinally stiffened cylinders. 
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Figure 9.- Predicted mode shapes for ring* and stringer-stiffened cylinder with if R - 0.15. 


